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Abstract: We compute the hadronic production of top-antitop pairs in association with 
a Higgs boson at next-to-leading-order QCD, including the decay of the top and antitop 
quark into bottom quarks and leptons. Our computation is based on full leading and 
next-to-leading-order matrix elements for e'’“Ue/U~r'^bbH(j) and includes all non-resonant 
contributions, off-shell effects and interferences. Numerical results for the integrated cross 
section and several differential distributions are given for the LHC operating at 13 TeV using 
a fixed and a dynamical factorization and renormalization scale. The use of the dynamical 
instead of the fixed scale improves the perturbative stability in high-energy tails of most 
distributions, while the integrated cross section is hardly affected differing by only about 
one per cent and leading to almost the same K factor of about 1.17. 
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1 Introduction 

After the discovery of a new boson of a mass around 125 GeV by the CMS and ATLAS 
Collaborations [1, 2] at the LHC, the determination of its quantum numbers and couplings 
to other particles has become a high priority in particle physics. Results from the first run 
of the LHC strongly support the hypothesis that this particle is the Higgs boson predicted 
by the Standard Model (SM) of particle physics. In the Brout-Englert-Higgs symmetry 
breaking mechanism of the SM the Higgs boson is the key for understanding the origin 
of mass. Since in this framework the Higgs boson couples to fermions with a strength 
proportional to their mass via Yukawa interactions, its coupling to the heaviest quark, the 
top quark, is of particular interest. The main production mechanism of the Higgs boson in 
the SM is gluon fusion, gg^H. This process is sensitive to the top-quark Yukawa coupling, 
but possible heavy particles beyond the SM running in the loop bias its determination. On 
the other hand, the production of a SM Higgs boson in association with a top-quark pair 
allows a direct access of the top-quark Yukawa coupling already at tree level, disentangling 
it from possible beyond-SM contributions. 

Leading-order (LO) predictions for ttH production for stable Higgs boson and top 
quarks have been presented in Refs. [3-7]. Since LO predictions of QCD processes suf¬ 
fer from large perturbative uncertainties higher-order corrections have to be taken into 
account for adequate theoretical predictions. The cross section for ttH production at next- 
to-leading-order (NLO) QCD is known for more than 10 years [8-12]. Meanwhile, NLO 
QCD corrections have been matched to parton showers [13-15] and recently electroweak 
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corrections to ttH production have been computed [16-18]. NLO QCD corrections for 
the important background processes ttbb and ttjj production have been worked out in 
Refs. [19-22] and Refs. [23-25], respectively, and matched to parton showers in Refs. [26-28] 
and Ref. [29]. 

Several searches for the associated production of the Higgs boson with a top-quark pair 
for a variety of decay channels have been published by ATLAS [30-34] and CMS [35-40]. 
These searches are challenging due to a large background from ttbb and ttjj production, and 
so far no evidence of a ttH signal over background has been found. The current ratio of the 
measured ttH signal cross section to the SM expectation quoted by ATLAS is = 1.5 ± 1.1 
[34] and by CMS ^ = 1.2^J5 [40] for a Higgs-boson mass of 125 GeV. 

In this article we present the calculation of the NLO QCD corrections to the hadronic 
production of a positron, a muon, missing energy, two b jets and a SM Higgs boson, which we 
assume to be stable (pp —at the 13TeV LHC, which includes the resonant 
production of a top-antitop-quark pair in association with a Higgs boson with a subsequent 
leptonic decay of the top and the antitop quark. Our calculation includes all NLO QCD 
correction effects in ttH production and top decays and also takes into account all off-shell, 
non-resonant and interference effects of the top quarks. We consider the fixed renormal¬ 
ization and factorization scale used in Ref. [9] and alternatively the dynamical scale choice 
from Ref. [13] and investigate their quality in reducing the dependence of the integrated 
cross section as well as differential distributions on the factorization and renormalization 
scales. The phase-space integration is performed with a newly implemented in-house multi¬ 
channel Monte Carlo program, using phase-space mappings similar to Ref. [41]. The Monte 
Carlo implements the dipole subtraction method [42-45] for the computation of the real 
corrections and is linked to the matrix-element generator Recola [46] for the computation 
of the LO and NLO matrix elements as well as colour- and spin-correlated squared matrix 
elements needed for the evaluation of subtraction terms. 

The calculation follows in many respects the one of pp —in Ref. [47]. 
In particular, since the additional Higgs boson in the final state has no colour charge, 
the contributing Catani-Seymour dipoles [42, 44] are the same. For the hadronic process 
pp —>■ e^r'e/.i~f^/ibbH we find a qualitative similar behaviour of NLO corrections to the 
integrated cross section as well as for differential distributions of the same observables. 
Where appropriate we compare results and point out differences. Moreover, we compare 
our results for the total cross sections to those of existing calculations for on-shell Higgs 
boson and top quarks [9, 13]. 

The paper is organised as follows. In Section 2 we discuss details of the calculation 
such as contributing subprocesses and Feynman diagrams and explain some technical as¬ 
pects of the real (Section 2.1) and virtual (Section 2.2) corrections. We present numerical 
results for the LHC operating at y/s = 13TeV in Section 3. In Section 3.1 we list and ex¬ 
plain the input parameters, jet definition, cuts and scale choices for our calculation, while 
results for integrated cross sections and a discussion of the scale dependence are provided 
in Section 3.2. In Section 3.4 we display and discuss several differential distributions for a 
fixed and a dynamical scale choice. We have performed several checks of our calculation 
which we present in Section 4. Our conclusions are given in Section 5. 
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Figure 1: Representative tree-level Feynman diagrams with (a)-(c) two (top), (d)-(f) one 
(middle) and (g)-(i) no top-quark resonances (bottom). 


2 Details of the calculation 


We compute the QCD corrections to the full hadronic process 


pp —)• r'^bbH. 


( 2 . 1 ) 


We consider the tree-level amplitude at 0(asa®^^) including all resonant, non-resonant, 
and off-shell effects of the top quarks and all interferences. Neglecting flavour mixing as 
well as contributions from the suppressed bottom-quark parton densities and counting u, d, 
c and s quarks separately, we distinguish 5 partonic channels for the LO hadronic process: 
the gluon-induced process gg —)• and four processes from qq 

by substituting different quark flavours (q = u, d, c,s). Throughout this paper we consider 
the bottom quark massless, implying no contribution from tree diagrams involving the 
Higgs-bottom-quark coupling. The gg process involves 236 and the qq processes 98 tree 
diagrams each under these prerequisites. In Figure 1 we show sample diagrams grouped by 
the number of top-quark resonances. 
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Figure 2: Representative hexagon and heptagon one-loop Feynman diagrams with two 
top-quark resonances. 

2.1 Real corrections 

The real correction process pp —)• e"’“r'e^~i^^bbHj receives contributions from the 13 partonic 
subprocesses 

gg ^ e+r'e/.i"i^AibbHg, 

qq e+r'elt"i^/.bbHg, 
gq e+t'e//“r';ibbHg, 
gq e+r'eAi“^'/^bbHg, 

where the gg process involves 1578 tree diagrams and the gg, gg and gg processes, all related 
by crossing symmetry, 614 tree diagrams each. 

Gluon Bremsstrahlung in the real corrections gives rise to IR divergences by soft or 
collinear configurations, which cancel for the final state for infrared-safe observables upon 
combination with the virtual corrections. Singularities from collinear initial-state split¬ 
ting factorize and can be removed by MS redefinition of the parton distribution functions. 
We employ the Catani-Seymour subtraction formalism [42, 44] for the regularization and 
analytical cancellation of IR singularities. Both the amplitudes for the real-correction sub¬ 
processes as well as the colour and spin-correlated amplitudes of the subtraction terms have 
been calculated with Recola. 

2.2 Virtual corrections 

The partonic subprocesses for the virtual QCD corrections can be identified with those 
at LO. We compute the virtual corrections in the ’t Hooft-Feynman gauge, where the 
gg process involves 9074 loop diagrams and the gg processes 2404 loop diagrams each. 
The most complicated one-loop diagrams are heptagons (Sample diagrams are displayed in 
Figure 2). 

The resonant top quarks, Z bosons and W bosons are treated in the complex-mass 
scheme [48-50], where the masses of unstable particles are consistently treated as complex 
quantities leading in particular to a complex weak mixing angle, 

/^w = “ i-^wrw, = M| — iMzTz, cos 9^ = (2.3) 

For the renormalization we use the on-shell renormalization scheme as described in Ref. [49] 
for the complex-mass scheme. 
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For the computation of the matrix elements for the virtual corrections we employ 
Recola [46] in dimensional regularisation, which integrates the COLLIER [51, 52] library 
for the numerical evaluation of one-loop scalar [53-56] and tensor integrals [57-59]. We com¬ 
pared our results for the virtual NLO contribution to the squared amplitude, 2ReA4QA4i, 
for many phase-space points with MadGraph5_aMC@NLO [60] (see Section 4 for details). 


3 Numerical Results 

3.1 Input parameters, jet definition, cnts and scale choice 


We present results for integrated cross sections and differential distributions for the LHC 
operating at y/s = 13TeV. For the computation of the hadronic cross section we employ 
LHAPDF 6.05 with CTIONLO parton distributions at LO and NLO QCD. We use the value 
of the strong coupling constant Og as provided by LHAPDF based on a one-loop (two-loop) 
accuracy at LO (NLO) with N-p = 5 active flavours. In the renormalization of Og the top- 
quark loop in the gluon self-energy is subtracted at zero momentum. The running of Og in 
this scheme is generated by contributions from light-quark and gluon loops only. For the 
fixed renormalization and factorization scale = 236 GeV, dehned below in (3.10), we 
find 


as(hfix) = 0.103237... . 


(3.1) 


We neglect contributions from the suppressed bottom-quark parton density. 

The electromagnetic coupling a is derived from the Fermi constant in the scheme 

[61], 

t ^ -- (3 2) 


a = —G,M? 


vr 


w 


M| 


Gf, = 1.16637 X 10-^GeV. 


We compute the width of the top quark Ft for unstable W bosons and massless bottom 
quarks according to Ref. [62]. At NLO QGD it reads 


pNLO _ 
^ t — 




dyyw 


16 V 27 r 2 M^ Jo (1 - y/yf 7^ 
with 7w = Fw/Mw, y = (Mw/mt)^, Og = as(mt)|NLO 


20- 

Fo{y)-^F^{y) 


(3.3) 


Fo{y) = 2{l-y)\l + 2y) (3.4) 

Fi(y) = 2(1 - yfil + 2y) [tt^ + 2Li2(y) - 2Li2(l - y)] (3.5) 

+ 4y(l - y- 2y2) ln(?/) -h 2(1 - yf{5 + Ay) ln(l - y) 

-( 1 - 2/)(5 + 92 /- 62 / 2 ). 


For the top-quark width at LO we neglect the correction term —2asFi{y)/[2>'K). 

As input, we employ the following numerical values for the masses and widths: 

rrit = 173 GeV, = 1.472886 ... GeV, = 1.346449 ... GeV, 

= 91.1876 GeV, = 2.4952 GeV, 

= 80.385 GeV, F^® = 2.0850 GeV, 

Mh = 126 GeV, 
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which includes the measured value of the Higgs-boson mass with zero width since we assume 
it to be stable. We neglect the masses and widths of all other quarks and leptons. 

We convert the measured on-shell values (OS) for the masses and widths of the W and 
Z boson into pole values for the gauge bosons {V = W, Z) according to Ref. [63], 

Mv = MOS/y/f + (rOS/MOS)2 , Ty = + (rOS/MOS) 2 , (3.7) 

that enter the calculation. 

We use the anti-^T algorithm [64] for the jet reconstruction with a jet-resolution pa¬ 
rameter R = 0.4. The distance between two jets i and j in the rapidity-azimuthal plane is 
defined as 

Rij = \J{<t>i - <i>jY + iVi - VjV, (3-8) 

with the azimuthal angle and the rapidity y* = ^ In of jet i, where E is the energy 
and pz the component of momentum along the beam axis. Only final-state quarks and 
gluons with rapidity |?/| < 5 are clustered into infrared-safe jets. 

After recombination we impose standard selection cuts on transverse momenta and 
rapidities of charged leptons and b jets, missing transverse momentum and distance between 
b jets according to (3.8). We require two b jets and two charged leptons in the final state, 
with bottom quarks in jets leading to b jets, and 


b jets: 
charged lepton: 
missing transverse momentum: 

b-jet-b-jet distance: 


PT,b>25GeV, |?/b| < 2.5, 
PT,t>20GeV, \yi\< 2.5, 
PT,miss 20 GeV, 

Ai?bb > 0.4. 


(3.9) 


We have identified the renormalization scale with the factorization scale p = /tr = PF 
and have considered a fixed reference scale set to half the partonic threshold energy for ttH 
production according to Ref. [9]: 

PFx = yR = y-F = ^ (2mt -f mn) = 236 GeV. (3.10) 

Alternatively, we use a dynamical scale following Ref. [13] 

Pdyn = PR = PF= {mT,t ® with mx = Jm? + p\, (3.11) 


which corresponds to the geometric mean of the top-quark, antitop-quark and Higgs-boson 
transverse masses. For the comparison of the scale choices we compute the logarithmic 
scale average /idyn of the dynamical scale, defined as 


In /idyn 


J ln(/idyn)do- 

/da 


(3.12) 


The scale uncertainty of the LO and NLO cross section is determined by variation of 
the renormalization and factorization scales pn and p-p around the central value pQ = ^ax 








Po 

ch. 

o-LO [fb] 

O'NLO [fb] 

K 

Mdyn 

gg 

1,5906(1)^1™ 

2.024(3)«« 

1.273(2) 


QQ 

0.67498(9)+24.i% 

0.495{1)^5™| 

0.733(2) 


(— ) 

SQ 


0.136(1)+S 



PP 

2.2656(1)+30.8% 

2.656(3)+°;^% 

1.172(1) 

/^fix 

gg 

1.5681(1)+“ 

2.011{3)(:S. 

1.282(2) 


QQ 

0.67199(9)+2+go 

0.495(1)+^+°! 

0.737(2) 


{—) 

SQ 


0.127(1)+?^ 



PP 

2.2401(1)+^iS 

2.633(3)+0;[;^“ 

1.176(1) 


Table 1: Composition of the integrated cross section for pp —)• e^r'e/^~i^/ihbH(j) at the 
LHC at y/s = ISTeV with both the dynamical (//dyn) and the fixed scale (^fix) as denoted 
in column one. In column two we list the partonic initial states, where q = u, d, c, s and 
Q = Qj Q- The third and fourth column give the integrated cross sections in fb for LO 
and NLO, resp. The upper and lower variations correspond to the envelope of seven scale 
pairs (hr/H o, Hf/ho) = (0.5,0.5), (0.5,1), (1,0.5), (1,1), (1, 2), (2,1), (2, 2). The last column 
provides the K factor with K = ciNLo/f^LO- 

and Ho = ^dyn for the fixed and dynamical scale choice, respectively.. While varying the 
renormalization scale in PDFs and matrix elements, the top-quark width remains fixed 
as computed at the top-quark mass according to (3.3). For the investigation of the scale 
dependence of the integrated cross section in Figure 3 we vary the scale h up and down 
by a factor of eight for the LO and NLO integrated cross section for the three cases: 1) 
Hr = t^F = H, 2) HR^ = H, t^F = Ho, 3) Hf = H, t^R = Mo- While we show all three cases for 
the dynamical scale, we show only the first case for the fixed scale choice in Figure 3. For all 
other results, i.e. those in Table 1 and Figures 5-7, the scale uncertainties are determined 
from factor-two variations as follows. We compute integrated and differential cross sections 
at seven scale pairs, {hr/H o, Hf/ho) = (0-5, 0.5), (0.5,1), (1,0.5), (1,1), (1,2), (2,1), (2, 2). 
The central value corresponds to {hr/H o, Hf/H o) = (1; 1) und the error band is constructed 
from the envelope of these seven calculations. 

3.2 Integrated cross section and scale dependence 

In Table 1 we present the integrated cross sections with fixed (3.10) and dynamical scale 
(3.11) at the LHC at y/s = 13TeV corresponding to the input parameters (3.1), (3.2) and 
(3.6) and the cuts as defined in (3.9). The results include only contributions of ©(osO^/^) 
for LO amplitudes and the corresponding 0(as) QCD corrections. We neglect possible 
contributions to qq processes of for LO amplitudes, which we determined to be 

about 2 per mille of the integrated cross section at LO for the setup described above. We 
also do not include partonic channels with incoming bottom quarks. At LO and using the 
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Figure 3: Scale dependence of the LO and 
NLO integrated cross section at the 13TeV 
LHC. The renormalization and factorization 
scales are varied around the central values 
of the fixed (/tq = dash-dotted lines) 

and dynamical scale (^uq = fJ-dyn, solid lines). 
For the dynamical scale the variation with 
/XR while keeping /xr = /Udyn fixed and vice 
versa is shown with dashed lines. 



r, [GeV] 

Figure 4: Zero-top-width extrapolation of 
the LO and NLO cross section at the LHC 
at y/s = 13TeV for fixed scale /xq = /^ax- 


dynamical scale (3.11), these contribute 4.869(1) x 10“^ fb, i.e. 0.2%, to the integrated cross 
section. Moreover, we calculated the integrated cross section at LO with finite bottom-quark 
masses, leading to a reduction of the cross section for our set of cuts by 0.03%. 

The use of the dynamical scale instead of the fixed scale increases the LO and NLO 
cross sections by only about 1%, and the K factor is 1.172 for the dynamical scale and 
1.176 for the fixed scale. The similar quality of both scale choices is also supported by the 
logarithmic scale average of the dynamical scale as defined in (3.12). With fidyn = 222.3 GeV 
it corresponds to a slight effective decrease of the fixed scale by only about 6%. 

Since integrated cross sections and NLO effects are very similar, the following con¬ 
siderations hold true for both, the dynamical and fixed scale: The major contributions 
to the cross section originate from the gluon-fusion process, with about 70 % at LO while 
increasing at NLO to 76 %. The contribution of the quark-antiquark annihilation drops 
from about 30% at LO to 19%. At NLO the gluon-(anti)quark induced real-radiation 
subprocesses contribute about 5 % to the integrated cross section. The inclusion of NLO 
QCD corrections reduces the scale dependence from 31 % to 5 %. 

We display the dependence of the integrated LO (blue) and NLO (red) cross sections 
on the values of the fixed and dynamical scale in Figure 3. Solid lines for the dynamical 
scale and dash-dotted lines for the fixed scale show the scale dependence for a simultaneous 
variation of the renormalization and factorization scales and dashed lines the individual 
















variation, where one of the scales is kept fix at the central value, for the dynamical scale only. 
While the largest scale variation is obtained when both scales are changed simultaneously, 
the smallest effect results if only the factorization scale is varied. The cross sections for the 
fixed and dynamical scale choices are uniformly shifted relative to each other by about 1 % 
as for the central scale /xq both for LO and NLO except for /x < /xo/2, where the fixed scale 
leads to a faster decrease of the cross section with /x as the dynamical scale. For the fixed 
and dynamical scale the maximum of the NLO cross section is near /x ~ justifying the 
use of both scale choices to be stable against scale variations. The K factor equals one at 
the slightly lower scale of about ^ ~ 0.7/xo 

3.3 Limit of on-shell top quarks 

To determine the effects of non-resonant and off-shell top-quark contributions on the in¬ 
tegrated cross section we perform a numerical extrapolation to the zero-top-width limit, 
Tt —>■ 0. To this end we plot 

^“/NLO(r.) = a“/"“(r.) j' ( 3 . 13 ) 

in the range 0 < Ft < where jg the top-quark width at LO and NLO, 

resp., and extrapolate linearly to Ft ^ 0, using a linear regression based on the computed 
LO and NLO integrated cross sections, as shown in Figure 4. The factor re¬ 

stores the physical top-decay branching fraction. Finite-top-width effects can be extracted 
by comparing the results for —)• 0) to (jLO/NLO^pLO/NLO^ LHC at 

^/s = 13TeV for fixed scale /xq = /Ufix finite-top-width effects shift the LO and NLO cross 
section by —0.07 ± 0.01 % and —0.14 ± 0.22 %, respectively, which are within the expected 
order of Ft/mt- The strong suppression of finite-top-width effects is related to the require¬ 
ment of a final state with two hard b jets. Finite-width effects are much more sizable in 
calculations where phase-space regions allowing for associated single-top plus W-boson pro¬ 
duction are included [65]. Such calculations require massive bottom quarks to regularize 
collinear singularities. 

3.4 Differential distributions 

In this section we present various differential distributions with two plots for each observ¬ 
able: The upper plot showing the LO (blue, dashed) and NLO (red, solid) predictions with 
uncertainty bands from the envelope of scale variations by seven pairs, /^^o) = 

(0.5,0.5), (0.5,1), (1, 0.5), (1,1), (1, 2), (2,1), (2, 2). The lower plot displays the LO (blue) 
and NLO (red) predictions normalized to the LO results at the central scale, i.e. TFlo = 
do'Lo(f^)/d(TLo(/^o) and LCnlo = dcTNLo(/^)/d(TLo(/^o)- Thus, the central red curve corre¬ 
sponds to the usual NLO correction factor {K factor), defined as K = o'NLo(/^o)/'7Lo(f^o)- 
The blue band shows the relative scale uncertainty of the LO differential cross section. 

Most of the displayed differential distributions were obtained using the dynamical scale 
/Xdyn, except for Figure 5 which illustrates the effect of the fixed-scale choice on transverse- 
momentum distributions. 
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PT,e+ [GeV] 



PT,bi [GeV] 


Figure 5: Transverse-momentum distributions at the LHC at ^/s = 13TeV for fixed scale 
Ho = (a) for the positron (left) and (b) for the harder b jet (right). The lower panels 

show the K factor. 


Where appropriate we compare NLO effects in differential distributions to those of the 
related process pp —)• e+z/e/U^i^/ibb presented in Ref. [47]. In general we find similar NLO 
effects for most of the distributions, but being often more distinct for pp —>■ e"*“r'e/^~f^^bb in 
Ref. [47]. 

Figures 5a and 5b display the transverse-momentum distributions of the positron and 
the harder b jet, resp., for the fixed scale ho = /^fix- The K factor of the transverse- 
momentum of the positron drops by about 50 % within the plotted range. In the high-px 
tail the NLO predictions move outside the LO band with a scale variation of almost the 
same size as the LO one. The K factor of the transverse-momentum of the harder b jet 
exhibits the same tendency, but not as drastic as for the positron. It moves outside the LO 
band for px < 60 GeV with larger scale variation as at the average px of around 90 GeV. 

In Figure 6 we collect several transverse momentum distributions obtained using the 
dynamical scale po = Mdyn- Figures 6a and 6c show the transverse-momentum distributions 
of the positron and the harder b jet, resp., to compare with the fixed-scale distributions 
in Figures 5a and 5b described above: They show clearly that the dynamical-scale choice 
improves the perturbative stability. The K factor changes only slightly (within 20 %) over 
the displayed range, and the NLO band lies within the LO band. The residual scale variation 
is at the level of 10% at NLO. 

In Figure 6b the distribution of missing transverse momentum, defined as px,miss = 
|pT,i/e + Pt,£>^|) is shown. The K factor rises for px,miss ^ 100GeV up to about 1.5. 

Figures 6c and 6d display the distribution of the transverse momentum of the harder 
and softer b jet, resp. While the K factor for the harder b jet exhibits a minimum at the 
maximum of the distribution and slightly rises towards its tail, the K factor of the softer 
b jet decreases by about 30 % in the plotted range. 

The distribution of the transverse momentum of the b-jet pair in Fig. 14 of Ref. [47] 
exhibits a strong suppression of the tt cross section at LO above Px bb ~ GeV. This 
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Figure 6: Transverse-momentum distributions at the LHC at y/s = 13TeV for dynamical 
scale /io = fJ-dyn- (a) for the positron (upper left), (b) for missing energy (upper right), 
(c) for the harder b jet (middle left), (d) for the softer b jet (middle right), (e) for the 
b-jet pair (lower left) and (f) for the Higgs boson (lower right). The lower panels show the 
K factor. 
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Figure 7: Differential distributions at the LHC at ^/s = 13TeV for dynamical scale /tq = 
fidyn- invariant mass of (a) the ttH system (upper left) and (b) the b-jet pair (upper 
right), (c) the cosine of the angle between the positron and the muon (middle left), (d) the 
azimuthal angle between the positron and the muon in the transverse plane (middle right), 
the rapidity of (e) the top quark (lower left) and of (f) the Higgs boson (lower right). The 
lower panels show the K factor. 
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is due to the fact that in narrow-top-width approximation the b quarks are boosted via 
their parent top and antitop quark, which have opposite transverse momenta resulting in a 
suppression of a bb system with high px at LO. The lesser stringent kinematical constraints 
at NLO result in an enhancement of the cross section and thus a large K factor for high 
PTbb- production at hand a Higgs boson acquiring transverse momentum 

softens the kinematical constraint already at LO leading to smaller NLO corrections for 
high Pxbb! which can be seen in Figure 6e. Furthermore, the K factor of the distribution 
resembles the K factor of the missing transverse momentum due to a kinematically similar 
configuration, but with a stronger increase to a value of 1.8 at px — 400 GeV. 

Figure 6f displays the transverse-momentum distribution of the Higgs boson. The av¬ 
erage px of the Higgs boson is around 70 GeV. The cross section decreases more moderately 
with px in the plotted range than for other transverse momentum distributions presented 
in Figure 6. 

In Figure 7 we present further differential distributions for other types of observables: 
the invariant mass of the ttH and bib 2 system in Figure 7a and Figure 7b, resp., the 
cosine of the angle between the two charged leptons in Figure 7c, the azimuthal angle 
in the transverse plane between them in Figure 7d, and the rapidity of the top quark 
and the Higgs boson in Figure 7e and Figure 7f, resp. Large NLO corrections appear in 
the Mttjj distribution below the ttH threshold. These arise dominantly from real gluon 
bremsstrahlung contributions, where the emitted gluon is not recombined with the bottom 
quarks and thus does not contribute to M^m. The distribution of the azimuthal angle 
in the transverse plane between the two charged leptons exhibits sizeable NLO effects for 
small angles similarly as the distribution in the cosine of the angle between the two charged 
leptons. The K factor varies by 40% for these distributions. The rapidity distribution 
of the Higgs boson features about the typical NLO effects in the central detector region, 
which disappear going into forward or backward direction. The K factor of the top-quark 
rapidity distribution on the other hand is almost flat at the level of the NLO effects of the 
integrated cross section. 

4 Checks 

We have performed several comparisons and consistency checks of our calculation. We 
have reproduced the LO hadronic cross sections with MadGraph5_aMC@NLO [60] using 
the fixed scale choice. We also compared the virtual NLO contribution to the squared 
amplitude, 2Re A4qA4i, for the gg and uu subprocesses computed by Recola [46] for many 
phase-space points with MadGraph5_aMC@NLO. In Figure 8a we plot the cumulative 
fraction of events with a relative difference larger than A between 2Re A4gA4i obtained by 
Recola and MadGraph5_aMC@NLO. The median of agreement is about 10“^ for the gg 
subprocess and roughly 5 x 10“® for the uu subprocess. The agreement is worse than 10“^ 
for less than 0.3 % of gg- and less than 0.04 % of uu-subprocess events. For the gg process we 
checked the accuracy of virtual matrix elements in satisfying the Ward identity, by replacing 
the polarization vector of one of the initial-state gluons in the one-loop amplitude A4i by 
its momentum normalized to its energy: Cg ^ pg/p^. In Figure 8b we plot the cumulative 


- 13 - 




Figure 8: Checks of the virtual contributions for fixed scale (a) agree¬ 

ment between Recola and MadGraph5_aMC@NLO for the partonic subprocesses uu —)• 
e'’“r'e/r~^'^bbH and gg —)• v (left), (b) accuracy in satisfying the Ward identity 

for gg —>■ e'’'r'e/r~^'^bbH (right), both showing the cumulative fraction of events with a 
relative difference/accuracy larger than A. 


fraction of events with 2Re AlQ(eg)Adi(eg ^ pg/pg)/2Re Adg(eg)Adi(eg) > A for virtual 
events obtained by Recola with a median of about 10“®. 

With the Monte Carlo code we developed for the calculation of the process pp —)• 
e'’'z/e/r~i^/ibbH, which employs Recola for the evaluation of matrix elements, we reproduced 
the results of Ref. [47] for the closely related process pp —>■ for the LHC at y/s = 

8TeV using a dynamical scale. Since the Catani-Seymour dipoles [42, 44] are the same for 
the process without the Higgs boson, this comparison allows us to test the implementation 
of the subtraction formalism. As we use the same renormalization procedure as in Ref. [47], 
this is also verified via this comparison. In Ref. [47] the one-loop matrix elements and the I- 
operator are computed in double-pole approximation for the two W-boson resonances using 
physical (i.e. real) W- and Z-boson masses and Lw = Fz = 0. Moreover, the Mh —)• oo 
limit is adopted, i.e. closed fermion loops involving top quarks coupled to Higgs bosons are 
neglected. In contrast, the new code computes the full one-loop matrix elements and the 
/-operator without any approximation. The largest difference between both calculations 
results from the use of the double-pole approximation for the virtual corrections and is of 
order ~ 0.1 %. We could reproduce the results of Ref. [47] for the integrated 

cross section at NLO within (2—3)(T, which is at the level of 0.5 %, thus confirming both 
calculations basically within the accuracy of the numerical integration. 

We compared our predictions for the NLO total hadronic cross section with compu¬ 
tations of ttH production without decays of the top quarks at NLO for fixed [9] and dy¬ 
namical scale [13]. To this end we performed NLO computations using the references’ 
parameters and PDF mappings and a minimal set of cuts on the decay products of the 
top and antitop in our process. Since our calculation includes background processes to 
ttH production, a minimal set of cuts has to be maintained to ensure infrared safety: to 
avoid large contributions from possible collinear events with bottom quarks in forward di¬ 
rection (as induced by diagrams shown in Figure Ig or Figure Ih) we kept a small b-jet 
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transverse-momentum cut (pT,b > 2 GeV) and a small b-jet distance cut (Ai2bb > 0.01) to 
avoid singular events from gluon splitting (g —)• bb) as in Figure li. We multiply appropri¬ 
ate branching ratios to the results of Refs. [9, 13] in the narrow-top-width approximation 
(NtWA), i.e. f ddNtWA = <7ttHBRt->-iBRt_^j, and apply the NLO matching prescription of 
Section 2.1.2 in Ref. [47] to our results. Thus, we find agreement of our NLO predictions 
with Ref. [9] and Ref. [13] within 1%, which is of the expected order of Ft/mt, since our 
calculation includes off-shell and non-resonant top-quark effects. 


5 Conclusions 

In this article we have presented the calculation of the next-to-leading-order QCD correc¬ 
tions to off-shell top-antitop production in association with a Higgs boson with leptonic 
decay of the top quarks at the LHC, including all resonant, non-resonant, and off-shell 
effects of top quarks as well as all interferences. For the computation of leading- and next- 
to-leading-order matrix elements we have utilised the recursive matrix-element generator 
Recola linked to the one-loop integral library COLLIER. The phase-space integration 
has been performed with an in-house multi-channel Monte-Carlo code that implements the 
dipole subtraction formalism. 

We provided integrated cross sections and several differential distributions for a 13 TeV 
LHC using a fixed and a dynamical scale that have both been used in the literature for the 
computation of NLO QCD corrections of ttH production with a stable Higgs boson and 
stable top quarks. We find almost the same integrated cross sections and scale dependence 
for both scale choices at leading order as well as next-to-leading order QCD, with a similar 
K factor of 1.172 and 1.176 for the dynamical and fixed scale choice, resp. However, the 
use of the dynamical scale instead of the fixed scale improves the perturbative stability in 
high-energy tails of most distributions, especially those of transverse momenta. Using the 
dynamical scale, we find K factors in the range 1.0—1.4 and residual scale uncertainties are 
the level or 10 % for distributions. For the integrated cross section an extrapolation to the 
zero-top-width limit has been performed, indicating non-resonant and off-shell top-quark 
effects below one per cent. While this effect is small, our calculation is also the first one 
to include NLO correction effects in the top-antitop-Higgs-boson production and the top 
decay processes. 

Besides its phenomenological relevance, this calculation demonstrates the power of the 
tools Recola and Collier in performing complicated NLO calculations. 
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